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Abstract 
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are used to obtain spectroscopic factors and asymptotic normalization coefficients, and they can 
serve as an input for low-energy reaction calculations. 
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I. INTRODUCTION 



One-nucleon spectroscopic overlaps, defined as the expectation value of the nucleon re- 
moval operator between states of nuclei differing by one particle, have been used in the 
treatment of processes in which one particle is added to or removed from the system 
In particular, these functions provide the interface between nuclear structure and direct re- 
action theories Traditionally, direct reaction theories of processes such as one-nucleon 
transfer, capture, knockout, or electron-induced proton knockout assume that the state of 
motion of only one particle is changed as a result of the projectile-target interaction, while 
the other particles are left essentially undisturbed. Under such a drastic approximation, the 
structure input needed to analyze experimental reaction data consists of single-particle-like 
functions of the involved nucleon inside the interacting nuclei, the one-nucleon overlaps. 

For decades, direct reactions have been the main tool for extracting spectroscopic factors 
(SFs) defined as norms of spectroscopic overlaps. Given the single-particle-like nature of 
these overlaps, SFs are often associated with single-particle-level occupancies and compared 
mainly with shell-model values. In practical applications, "experimental" SFs are extracted 
from experimental cross sections in a mo del- dependent way with structure and reaction 
uncertainties entangled. This dependence was illustrated, for example, in 0, Isl, where 
SFs from transfer and electron-scattering data were shown to differ unless extracted in a 
consistent way. Strictly speaking, however, the overlaps and the SFs are merely theoretical 
concepts that cannot be measured The overlaps should not be interpreted as the measurable 
probability amplitudes of finding a nucleon at some distance from the rest of the nucleus Q. 
Nevertheless, until more advanced reaction models become available, spectroscopic overlaps 
and factors need to be calculated theoretically for a chosen Hamiltonian to enable the analysis 
of direct reaction data. These concerns are less severe for the long-range scaling factors of 
overlaps, the asymptotic normalization coefficients (ANCs), which enter the cross section of 
peripheral direct reactions, as these coefficients remain invariant under finite-range unitary 



transformations of the Hamiltonian [10[. The spectroscopic overlaps and factors are still 
useful to assess the relative importance of clustering patterns in a given nucleus. 

Typically, spectroscopic overlaps are approximated by shell-model or mean-field single- 
particle functions, or they are taken as solutions of single-particle (Woods-Saxon (WS)) 
potential wells with some commonly accepted but otherwise arbitrary potential parame- 
ters [2I, 3, 11|. Of a more recent vintage are some calculations attempting to account for 
contributions from missing model space, for example results obtained within the correlated 
basis function theory l^j or overlaps generated as solutions of an inhomogeneous equation 
with a shell-model source term [l3|. Realistic calculations of one-nucleon overlaps based on 
modern nuclear interactions are scarce and limited to light nuclei as they are complicated 
even in the simplest cases: in the s-shell, they were computed in the hyperspherical harmon- 
ics (HH) method [13, fisf : some overlaps and SFs for p-shell nuclei were calculated within 
the no-core shell modeljlG. 17 | and the coupled-cluster method [is!]. Of particular relevance 
for the present work are previous variational Monte Carlo (VMC) overlap calculations for 
nuclei up to A = 10 (unpublished, available from [l9|]) and a more recent VMC-based cal- 
culation of ANCs in light nuclei |20|. VMC overlaps have been used in some analyses of 
hadronic 21- and leptonic [26| experiments. 

This work provides a systematic study of one-nucleon spectroscopic overlaps, SFs, and 
ANCs in light nuclei calculated within the Green's function Monte Carlo (GFMC) method 
based on realistic two- and three-nucleon interactions. The GFMC method is designed to 
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project the exact solutions out of trial wave functions by propagating them in imaginary 
time. The GFMC propagations are initiated by VMC wave functions. Over the years, the 
VMC/GFMC method has been found to accurately describe the structure and some reaction 



aspects of light nuclei |27H33l|. Given our experience, GFMC is expected to improve VMC 
overlaps, especially in the p-shell. In this paper, we present GFMC overlap results for nuclei 
up to A = 7; heavier nuclei up to A < 10 will be the subject of a forthcoming paper. 

This paper is organized as follows. In Sec. [ITl the theoretical framework is established, 
the VMC and GFMC methods are briefly reviewed, and all overlap-related quantities are 
defined. Sec. IIIII contains technical and computational details and an error analysis of the 
GFMC overlap calculations. Results are presented and discussed in Sec. [IVl final remarks 
and conclusions are given in Sec. |Vl 



II. QUANTUM MONTE CARLO 

To calculate spectroscopic overlaps, we first construct wave functions \E'( J'^, T, Tj) for 
nuclei of interest as solutions of the nonrelativistic many-body Schrodinger equation: 

Hm{J\T,T,) = Em{J\T,T,), (1) 

with J'^, T, and Tz denoting the total spin-parity, isospin, and isospin projection, respec- 
tively. The y4-body Hamiltonians used in this work have the form 

AAA 

H = Y,Ki + Y,Vv+ E ^^^^^ (2) 

i i<j i<j<k 

where Ki is the nonrelativistic kinetic energy and Vij and Vijk are two-nucleon (2A^) and 
three-nucleon (3A^) potentials, respectively. Most of the results presented in this paper have 
been obtained for a combination of Argonne fig (AV18) [l^l 2A^ and Illinois-7 (IL7) [35] 3A^ 
realistic interactions. For testing and benchmarking purposes, some calculations involved 
AV18 supplemented by the older Urbana IX (UIX) [36] 3A^ force. The kinetic energy of 
the center of mass is zero because the wave functions used in this work are translationally 
invariant. 



A. Variational Monte Carlo 



The wave functions are constructed in two steps. First, a trial VMC approximation of 
is written and optimized by minimizing the energy expectation value as computed by 

Our latest variational functions have the form 138L 



Metropolis Monte Carlo integration [37 
Eq. 3.13] 



1^. 



l + Ui 



A 



3N 
ijk 



(3) 



The Jastrow wave function, is a fully antisymmetric state having the desired J^,T,Tz 
quantum numbers of the state of interest; the Uij and UfJ^ are two- and three-nucleon short- 
range correlation operators induced by dominant parts of 2N and 3A^ forces, respectively, 
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while iS is a symmetrization operator that restores the antisymmetry violated by the non- 
commuting character of Uij and UfJ^. At long range, appropriate boundary conditions are 
imposed on \E'y (39| . 

The Jastrow wave function for s-shell nuclei consists of a simple product of central pair 
and triplet correlation functions and an antisymmetrized spin-isospin state. For p-shell 
nuclei, we use two types of \Ef j: a shell-model-like trial function, which we call Type I, and 
a clusterized version denoted Type II. These two kinds of wave functions differ primarily in 
their treatment of correlations between p-shell particles and the long-range asymptotics. 

In Type I trial functions, the single-particle radial functions of p-shell particles are ex- 
ponentially decaying solutions of WS potential wells centered at the s-shell a core with 
potential parameters subject to a variational search. Pair correlations between s- and p- 
shell particles and among p-shell particles themselves are chosen to allow clusterization in 
the p-shell. Being less restrictive on long-range cluster decomposition compared with Type II 
functions. Type I functions are more widely applicable. Details of Type I functions for nuclei 



with A = 6, 7 can be found in 40 



Type II functions that may be used for nuclei with low-lying cluster break-up thresholds 
impose a strict cluster-cluster asymptotic decomposition. Examples are the ground states 
of ^Li and ^Li that may asymptotically decouple into a + d and a + t, respectively. In 
Type II trial functions, the p-shell single-particle radial functions are again solutions of p- 
wave differential equations with potentials containing WS and Coulomb terms, but with an 
added Lagrange multiplier that turns on at long range. The role of the Lagrange multiplier 
is to impose the cluster boundary condition: 

^y(r ^ oo) oc tpa^xW.^^i+i/2{2kr)/r, (4) 

where x = d or t for ^Li and ''Li, respectively, and r and / are the relative a-x distance 
and orbital angular momentum. The Sommerfeld parameter t] and the wave number k are 
related to the a-x separation energy in a given nucleus, which is set to its experimental value. 



Here, W is the Whittaker function defined in Eq. (13.1.33) of [4l|. The correlations between 



p-shell particles in ^Li and ''Li are adopted from exact deuteron and variational triton wave 



functions, respectively. More details on Type II trial functions are given in [32|, |42, |43[ . 

For either type of trial function, a diagonalization is carried out to find the optimal values 
of mixing parameters for states of different spatial symmetries among p-shell particles [27| . 



B. Green's function Monte Carlo 

Being variational solutions of Eq. ([1]), the VMC wave functions can be improved fur- 
ther by action of the operator limT-_!.oo exp [— {H' — Eq) r], which, for a given set of quantum 
numbers, projects out the exact lowest-energy state ^'o of a possibly simplified version H' 
of the desired Hamiltonian H . In practice, the operator is applied in small increments of 
the imaginary time r up to some finite value to yield a propagated wave function \['(r): 

^(r) = e-(^'-^«)"^v (5) 

Obviously, \E'(r = 0) = \I/y and \E'(r — )• oo) = \E'o. In practice, Eq. (jS]) is turned into an 
integral equation involving Green's functions with integrations performed by Monte Carlo 
methods, hence the name Green's function Monte Carlo. The energy Eq is an approximate 
guess for the true eigenenergy corresponding to ^'o- 
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For a given state of the nucleus, quantities of interest are evaluated in terms of a "mixed" 
estimate between and \l/(r): 



(^(r)|Q|^y) 

The desired expectation values would have ^I^(t) on both sides; by writing ^^(t) = \E't+5^(''") 
and neglecting terms of order [5\l/(r)]^, we obtain the approximate expression 

(0(r)) = ^^^^^'^j^^)^^^^ ^ (0(r)),, + [(0(r))M - (O)v^], (7) 

where {0)v is the variational expectation value. The expectation value of H' is an exception 



to Eq. ([7]) since for this operator the mixed estimate already gives the correct value [27 
For off-diagonal matrix elements between two different A-body nuclear states, Eq. can 
be generalized to a form involving two mixed estimates (32I . 33l . 



We note that the GFMC propagator involves a simplified Hamiltonian H' = AV8' + H^j^ 
based on a reprojection AV8' of the full AV18 two-body potential j4o|. Therefore, the GFMC 
wave functions are really eigenstates of H' rather than H, which, in general, should not be 
a problem given that H' is a good approximation of H. In this work, the three-nucleon 
part H'^j^ is either UIX or IL7, with the strength of their central repulsive parts adjusted so 
that {H') ^ (H). Although energies can be corrected perturbatively by adding {H — H') to 
{H'), this kind of correction is not possible for other observables. In what follows, we will 
be referring to the desired Hamiltonian H remembering that, in fact, we use H' in GFMC 
propagations. 



For a detailed description of the nuclear GFMC method, see the review j27j and references 
therein. 



C. Spectroscopic overlaps and factors, and asymptotic normalization coefficients 



Let us consider a decomposition of an A-body nucleus, the parent, into the core C 
and a valence particle v (a proton or a neutron), namely A ^ C + v. The terms "core" 
and "valence" are used in a figurative sense considering that the parent wave functions 
are explicitly antisymmetrized states of indistinguishable particles. Then, a one-nucleon 
spectroscopic overlap is defined as a projection of a parent state onto an antisymmetrized 
core-valence product form: 



(8) 



A, 



Cv 



S{r - rcv 



Ja 



'Cv 



where CJ'*^''''*.. 



is an isospin Clebsch-Gordan coefficient; a = {A, J^, T4, ^j, 7 = 
{C, Jq,Tc,Tz^c}, and u = {v,l, s, j,t,tz} are composite indexes for the parent, the core, 
and the valence particle, respectively; ^ comprise all degrees of freedom in a given nucleus; 
fcv is the vector extending from the core's center of mass to the valence particle; and Acv 
is the core- valence antisymmetrizer: Acv = ~^(~1)^-Pcv with Pcv being the 

permutation operators exchanging the valence particle with those inside the core and p 
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being the permutation parity ±1. The symbols A and C = A — 1 also denote the nuclear 
mass numbers, and the subscripts A and C may be dropped when convenient. In Eq. ([8]), 
the parent and core states are assumed to be normalized to unity, and the integral is done 
over all 3A dimensions. The valence angle-spin-isospin function is defined as 



(9) 



where is a spherical harmonic, Xs{cry) and Xt,tA'^v) denote the valence spin and isospin 
states with s = t = 1/2, and gives the valence isospin projection ±1/2. Taking into 
account the antisymmetry of the parent state, one can turn the integral in Eq. into a 
simpler form: 



R{a,-f,jy;r) 



A([^{^)®yiu)] 



Ja 



Sir - rcv) 



'Cv 



(10) 



By definition ([8]) and by using translationally invariant wave functions the overlap func- 
tions R are translationally invariant. For a given parent-core combination, different angular 
momentum channels will be denoted L. 



The theoretical SF is then defined as the norm of the overlap: 



S'(a,7, z/)= / |_R(a, 7, 1/; r)|^ dr. 



(11) 



Among other sum rules, all possible (proton, neutron) SFs for a given state of the parent 
nucleus add up to the parent's number of nucleons (protons, neutrons) |l[. Our definition 
of SFs is consistent with those of some other works, for example [l3| but differs from oth- 
ii by inclusion of the isospin Clebsch-Gordan coefficient in Eq. ([8]) into overlaps and 



ers 



consequently into SFs (these works define SFs as norms of 0{r) in Eq. ([8])). 

To calculate the overlap functions R within GFMC, one can, neglecting terms of or- 
der [\l/(Q;;r) — '^vi'^)]'^ and [\l/(7;r) — \l/y(7)]^, derive the following expression similar to 
Eq. (19) in Q: 



i?(a,7,z/;r;r) 



[^(7;r)®3^(z/)] 



Ja 



S{r~rcv) 
^ 



{R{a, 7, z/; r; t))ma + {R{a, 7, r; t))mc - (^(«, 7, ^] r))v, 



where 



(/?(«, 7, z/;r;r))MA = 
(i?(a,7, z/;r;r))Mc = 

(i?(a,7, z/;r))v = 
with a self-normalization factor: 



[*y(7)®3^(^)]j. 



<5(r-rci,) 
71 



' Cv 



T 





(a;r)) 

&{r-rcv) 




(^(7;r)h 


&y(7)) 

5{r-rcv) 


^y(a)) 




Af 



{<i/v{a)\^v{a)) 



Af 



(12) 

(13) 
(14) 
(15) 

(16) 
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The overlap function {R)v in Eq. ( !T5|) is a pure variational estimate, to be called a VMC 
overlap, whereas the other two estimates, {R)ma Eq. (fT3|) and {R)mc Eq. (fT4l) . involve 
combinations of GFMC-propagated and VMC wave functions, and as such will be called 
A- and (74-l)-body mixed overlaps and referred to by the corresponding GFMC-propagated 
nucleus. For each estimate Eqs. (IT^ - (IT^ . one can define a SF similar to Eq. ( ITT]) . The 
GFMC-extrapolated overlap R in Eq. fll2p will be called a GFMC overlap. Obviously, at 
T = 0, R = {R)v = {R)ma = {R)mc- 

Overlap functions R satisfy a one-body Schrodinger equation with the appropriate source 
terms j2|. Asymptotically, for r — )■ oo, these source terms contain core- valence Coulomb 
interaction at most, and hence the long-range part of overlap functions for parent states 
below core- valence separation thresholds is proportional to a Whittaker function W^ri,i+i/2'- 



R{a, 7, z/; r) C{a, 7, u) 



iy_^,;+i/2(2fcr) 



(17) 



where r] = ZqZ^ {e^ /he) fic^/2B depends on proton numbers Zq and Z^, and the core- 
valence reduced mass n and separation energy B (positive for parent states below core- 
valence separation thresholds). The wave number k is defined as \/2^B/fi, and / is the 
orbital angular momentum from Eq. Note that W , /i, B, and k implicitly_carry channel 
labels a, 7, and v. The Whittaker function is defined in Eq. (13.1.33) of [41], and it has 
an (approximately) exponentially decaying tail. The proportionality constant C(a,7,z/) in 
Eq. ([T7D is the ANC. 

To provide a convenient parametrization of overlaps suitable for reaction calculations and 
to extract ANCs, we perform x^-fiis of the overlaps by the eigenstates of a single-particle-like 
Hamiltonian — [/i^/2/i]A + V{r) containing a WS plus spin-orbit (so) plus Coulomb (Coul) 
potential 



V{r) = H 



WS 



1 



41 -s 



1 + exp((r - Rws)/aws) 
K„ d 



r dr 



1 + exp((r - Rso)/aso) 



- /3exp(-(r/p)2) 



+ Vcoul 



+ 



(18) 



with 



Coul 



ZcZ^e^r 
ZcZ^e^ [3 



[r/Rcouif] /[2Rcoui] 



■ r > Rcoui 
: r < Rcoui 



(19) 



The central WS part of the potential includes a Gaussian "wine-bottle" term to provide an 
additional flexibility at short range. The depths Vws and Vso are in the same units, the factor 
4 in the spin-orbit part is approximately twice the square of the pion Compton wave length 
^/(m^rc) in fm, and / and s'are operators of the core- valence orbital angular momentum and 
of the valence spin, respectively, both in units of Ti. The Coulomb radius is Rcoui = 2 fm. 
The potential parameters are varied freely to provide the best fit under the constraint of 
the eigenenergy being equal to the desired value —B. This fitting procedure provides good 
overlap fits at short and medium distances; at large distances, the fits have the desired form 
shown in Eg . ([T7|) . These overlap parametrizations can be used in reaction codes such as 
Ptolemy 4^ or Fresco jiej, and they will be refered to as WS fits throughout the paper. 
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III. INNER WORKINGS OF OVERLAP CALCULATIONS AND ERROR ANAL- 
YSIS 



In this section, we elaborate on some technical aspects of the overlap calculations, show 
several detailed examples, and assess systematic errors of quantities being computed. The 
discussion of results is given in Sec. HVl 

The VMC overlap {R)v in Eq. fllSp is calculated on a random walk guided by the parent's 
|^y(a)p. On the same walk, we also evaluate the self- normalization factor from Eq. fll6l) : 
to extend the denominator of Af into the 3y4-dimensional space being sampled, \&y(7) is 
replaced by "^vil) x pi^cv), where p, normalized to unity as / p'^{rcv) dfcv=i, acts as a 
single-particle function of a virtual valence particle taken with respect to the core's center of 
mass. In principle, any function can be used for p as long as the product form |\E'y(7) xpp is a 
reasonably good approximation of the A-hodj sampling probability density in order to yield 
reasonably small statistical errors on the denominator of A/". In practice, we find it sufficient 
to use purely radial (s-wave) functions p{rcv) of either a Gaussian shape or one generated by 
a single-particle-like Hamiltonian containing a WS potential well with parameters adjusted 
to minimize the error on the self-normalization factor. The WS-generated single-particle 
functions have the advantage of better approximating the A-hodj sampling density at large 
core- valence distances, especially when the potential depth is set to approximately reproduce 
the (experimental or VMC) core-valence separation energy. In the present work, the self- 
normalization factors A/" were calculated (after accounting for auto-correlations between local 
samples) with an accuracy of the order of 0.1% or better. 

The two mixed overlaps in Eq. (fT3|) and Eq. (IT^ are calculated on random walks guided 
by GFMC propagations for the corresponding nucleus. For the (A-l)-body mixed over- 
lap, the GFMC sampling density only spans a 3(A-l)-dimensional subspace of the full 
3A-dimensional space of the parent nucleus; we draw the position of the valence nucleon 
from the same single-particle density p'^{rcv) used in the computation of the variational 
self-normalization factor A/". 

Generating new VMC or GFMC samples is computationally more expensive than eval- 
uating a sample's contribution to Eqs. (fT3|) - (fT5|) . Therefore, in order to improve the com- 
putational efficiency, each VMC or GFMC sample is used several times: in Eq. f|T3l) and 
Eq. fllSp . we consider all A cyclic particle permutations to rotate the valence particle over 
all possible positions within the parent nucleus, while in Eq. (1141) the position of the valence 
particle is drawn several (A) times from its single-particle distribution p"^ described above 
for each core's GFMC sample. 

We illustrate the method with the (^H | ^He) overlap for AV18+IL7 taken as a typical 
example of overlaps in the s-shell. First, the VMC overlap is computed by Eq. (IT^ . Then, 
the wave functions of the parent and the core are independently propagated and the two 
mixed and the GFMC overlaps are computed in each radial bin following Eqs. 0131) . fll4p . and 
f|T2|) . The corresponding SFs are plotted in Fig.[Tl In all figures, the quantities corresponding 
to the variational, the (^4-1)- and the A-hodj mixed, and to the GFMC overlaps are plotted 
as black, red, blue, and green, respectively. 



For s-shell nuclei, we do unconstrained propagations |40| and thus obtain essentially 
exact solutions for a given Hamiltonian H' in Eq. (E]). Because of the fermion sign problem, 
the statistical noise grows with time, as can be inferred from Fig. [H Fortunately, GFMC 
quickly eliminates excited-state impurities from VMC wave functions resulting in a rapid 
convergence of mixed and GFMC overlaps and SFs. For s-shell nuclei, GFMC propagations 
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FIG. 1. (Color online) Imaginary time evolution of SFs in the Si/2 channel of the (^H | ^He) overlap 
obtained for the AV18+IL7 potential. Only statistical errors are shown. Horizontal lines are SFs 
for VMC and time-averaged mixed and GFMC overlaps, as described in the text. 



are fully converged long before r = 0.5 MeV , at which point they are terminated to avoid 
the increasing statistical noise. 

In order to improve the statistical accuracy further, the [A-l)- and A-body mixed overlaps 
are separately averaged over time. These mixed time-averages are then combined with the 
VMC overlap to obtain the final time-averaged GFMC overlap from Eq. (fT2|) . In Fig. [1] 
and similar figures, the variational estimate, although computed at r = MeV~^, is plotted 
across the full range of times, while the other solid lines give SFs for time-averaged mixed 
and GFMC overlaps with averaging done over the range of times indicated by the horizontal 
extent of these lines; the dashed lines show statistical errors. We have checked that our 
time-averaged results have little dependence on the exact time interval being averaged over 
as long as this interval is safely within the GFMC converged region and does not include 
large r's for unconstrained propagations. 

In Fig. m the VMC and time-averaged mixed and GFMC (^H | "^He) overlaps are shown 
along with the asymptotic Whittaker function corresponding to the experimental separation 
energy (discussed below), a WS fit to the GFMC overlap, and the radial sampling probabil- 
ity. In general, the differences between s-shell VMC and time-averaged mixed and GFMC 
overlaps in a given overlap channel are very small because the starting VMC wave functions 
are already very good approximate solutions of Eq. ([T]); these small differences are reflected 
by small (~2%) differences between VMC and GFMC SFs as illustrated in Fig. [T] Statisti- 
cal errors on SFs are small because the dominant contribution to Eq. f llip comes from the 
volume region well covered by the Monte Carlo sampling probability, also shown in Fig. [2J 

Overlaps between p-shell nuclei are computed by using the same algorithm with some 
technical modifications. Compared with s-shell nuclei, the fermion sign problem becomes 
more severe in the p-shell, and we retreat to constrained path GFMC sampling jssf. Conse- 
quently, statistical errors are well under control, as illustrated in Fig. Elfor SFs in the P3/2 
channel of the (^He(O^) | ^Li(3/2~)) overlap. Given our experience indicating that mixed 
estimates of many observables tend to fluctuate more in the p-shell than in the s-shell, we 
chose to carry out propagations for longer times up to r = 3 MeV~^ to ensure full conver- 
gence in the p-shell. VMC, and time-averaged mixed and GFMC overlaps corresponding to 
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FIG. 2. (Color online) Linear (left) and logarithmic (right) plots of VMC and time-averaged mixed 
and GFMC Sx/2 (^H I "^He) overlaps obtained for the AV18+IL7 potential. Only statistical errors 
are shown. Also shown are a WS fit to the GFMC overlap and the asymptotic Whittaker function 
corresponding to the experimental separation energy in the right panel. Superimposed in the linear 
plot is the sampling probability (arbitrary scale). 



Fig. [3] are shown in Fig. |H The difference of about 8% between the VMC and GFMC SFs 
reflects the change in the shape of the underlying overlaps. 

At short range, the GFMC overlap is a result of fine cancellations between VMC and 
mixed overlaps in Eq. f lT2|) . At large core- valence distances, the GFMC overlaps tend to 
follow the A-hodj mixed overlaps, while the (y4-l)-body mixed overlaps usually stay close 
to the VMC ones, in agreement with our expectation that it is the long-range fall-off of 
the parent's wave function that primarily sets the tail of the overlap whereas the detailed 
structure of the core plays a less important role in this region although the core still needs 
to be described reasonably well. From this point on, mixed overlaps and estimates will not 
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FIG. 3. (Color online) Imaginary time evolution of SFs in the P3/2 channel of the 
('5He(0+) I ■^Li(3/2")) overlap obtained for the AV18+IL7 potential. The GFMC wave functions 
originated from Type I VMC wave functions. Only statistical errors are shown. Horizontal lines 
are SFs for VMC and time-averaged mixed and GFMC overlaps, as described in the text. 
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FIG. 4. (Color online) Linear (left) and logarithmic (right) plots of VMC and time-averaged mixed 
and GFMC P3/2 (^He(0+) | ■^Li(3/2-)) overlaps obtained for the AV18+IL7 potential. The GFMC 
wave functions originated from Type I VMC wave functions. Only statistical errors are shown. Also 
shown are a WS fit to the GFMC overlap, and the asymptotic Whittaker function corresponding 
to the experimental separation energy in the right panel. 



be shown, and the time-averaged GFMC overlap will be called a GFMC overlap. 

Because of the vanishing sampling probabilities at large distances (an example is shown 
in Fig. [2]), reliable sampling of overlap tails requires very large Monte Carlo sets; we prop- 
agated about 12 million samples in the s-shell and about 250,000 for p-shell nuclei. As a 
consequence, statistical errors are small. It is therefore important to assess systematic errors 
in addition to statistical errors. The most notable sources of systematic errors are: the dif- 
ference between the Hamiltonian H' in the GFMC propagator in Eq. ([5]) and the desired H 
in Eq. ([1]), possible errors due to the constrained path sampling in the p-shell, and a residual 
dependence of GFMC results on starting VMC wave functions. In addition, ANCs depend 
on the separation energy used in their determination. We now attempt to place limits on 
such systematic errors. 

For s-shell nuclei, we can compare our results with those obtained by the hyperspherical 



harmonics (HH) method [IJ, ll5|, |47|, |48[ . For AV18+UIX, the agreement between HH and 
GFMC is very good, as can be seen from Fig. [6] for overlaps and from Table |IV] for SFs 
in Sec. IIVI Testing the bias due to H' for p-shell overlaps is not possible since, to our 
knowledge, no p-shell overlaps have been published previously for the realistic interactions 
employed here. To test the dependence of GFMC SFs on VMC wave functions, we include 
in Table |V] results involving the ground state of ^Li obtained for several combinations of 
Type I and II (see Sec. Ill Ap trial wave functions. Although VMC SFs for different 
may differ by as much as 20%, GFMC reduces the spread to no more than 3%. Hence, we 
estimate systematic errors on GFMC SFs to be no more than 2-3%. 

Typically, extracting ANCs directly from overlaps in Eq. fIlOp is problematic for most 
many-body methods. First, these methods may not yield the correct exponential asymptotic 
form given in Eq. (1170 . as is the case of methods employing harmonic oscillator bases |16l.ll7 . 
In variational methods, including VMC, it may be hard to impose the correct asymptotics 
while preserving short-range properties. This problem is less prevalent in integral methods 
of computing ANCs in which the desired asymptotics is imposed by definition regardless of 
the actual asymptotic shape of the wave functions involved 13, 13]. Second, ANCs extracted 
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r [fm] r [fm] 



FIG. 5. (Color online) The ratio C{r) of VMC, GFMC, and fitted GFMC overlaps to asymptotic 
Whittaker functions corresponding to experimental core-valence separation energies. Left panel: 
(^H I ^He) overlaps from Fig. [SJ right panel: (^He(0"'") | ^Li(3/2~)) overlaps from Fig. [H 



from Eq. (ITTj) depend on the separation energy B. To determine ANCs from overlaps in 
a fully consistent manner, one should use the separation energy Bh obtained for a given 
Hamiltonian; on the other hand, for ANCs to be practically usable for reaction calculations, 
experimental values B^xp should be used. Ideally, Bh ~ B^xp, a condition that is often 



violated in literature and therefore relaxed without error analysis for ANCs |l3l . llTl . |20| . We 
now discuss these aspects applied to GFMC. 

It is gratifying to observe that the GFMC overlaps in Fig. [2] and Fig. H] transition nicely 
from the interior into the asymptotic region and that in the computationally safe region 
(r < 8 fm) they follow closely the right asymptotic trend represented by the Whittaker 
functions. To illustrate this in more detail for overlaps from Figs. [2] and HI we show in Fig. [5] 
the ratios 

C(r) = R{r)/ [l^_,,i+i/2(2A;r)/r] . (20) 

Ideally, C(r — )■ oo) should approach C from Eq. f[T7|) . the ANC. As can be seen in the 
figure, the GFMC curves may not flatten out perfectly at larger distances where the method 
being driven by H' in Eq. ([5]) may break down because of insufficient contributions of the 
asymptotic part of the configuration space to the energy. To correct this imperfection, the 
GFMC overlaps are extended to larger distances by means of WS fits (Eq. (fT8|) ). For s- 
shell overlaps, C(r) fiatten out at r~2.5-3 fm where C(r)/C~0.9, and they are almost fully 
converged at r~5 fm; for p-shell overlaps, C(r)/C~0.9 is realized at r~5 fm. Therefore, the 
overlaps are (almost) asymptotic at about 5 fm, which is the upper radial limit used in our 
fitting procedure. ANCs are then determined by applying Eq. (1201) to GFMC overlap fits 
at large distances. We have checked that including the region beyond 5 fm in the fitting 
procedure does not result in significant changes of ANCs. 

In GFMC, the binding energy is thought to be computed with systematic errors of about 



1-2% [271, 138|, 140]. Most GFMC overlaps in this work were obtained for the AV18+IL7 
potentials. For this potential, GFMC energies in Table [I] differ by less than 1% from ex- 
perimental values, in most cases the difference is less than 0.5%. These deviations, though 
being small fractions of total binding energies, may translate into bigger (fractional) devia- 
tions on core- valence separation energies, as can be seen from Tables [TTl and IIIII Even then, 
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however, the difference between experimental and GFMC separation energes is ^ 3%, with 
a noticeable exception being the (^He(0"^) | ^He(3/2~)) overlap involving a particle- unbound 
^He (for this overlap, the ANC is not a meaningful quantity). Since GFMC wave functions 
are really eigenstates of H' rather than H, it is equally important to observe in Tables HIl 
and mil that separation energies for these Hamiltonians are about the same {H' was ad- 
justed to get {H') i^H) as mentioned in Sec. Ill Bp . To estimate the sensitivity of ANCs 
on separation energies, we refitted the GFMC overlaps from Figs. [2] and H] with solutions of 
Eq. f|T8l) imposing separation energies corresponding to ±3% deviations from B^xp- Using 
Bexp from Tables [TT] and IIIIl we obtain the following ANCs for B = 0.97Bexp, B = B^^p, 
and B = LOSBexp: -6.13, -6.45, -6.78 for the (^H | ^He) overlap, and 3.33, 3.52, 3.71 for 
the (^He | ^Li) overlap. Thus, any uncertainty of about 3% or less on the separation energy 
translates into an ANC uncertainty of about 5-6% or less. Furthermore, the s-shell GFMC 
ANCs in Table IVII agree well with those obtained within the HH method except for the 
weak d-waves in A = 3 nuclei for which the two methods differ for unknown reasons (the 
HH ANCs assume the actual separation energies for a given Hamiltonian that are within a 
fraction of percent of the experimental values). To test the dependence of GFMC ANCs in 
the p-shell on starting VMC wave functions, we include in Table IVIII results involving the 
ground state of ''Li obtained for several combinations of Type I/II (see Sec. Ill AO trial wave 
functions; depending on the starting VMC wave function, GFMC ANCs differ by at most 
7% in a given overlap channel. 

Based on these arguments, we chose to determine the GFMC ANCs assuming experi- 
mental separation energies. We estimate the systematic errors on our ANCs to be ^ 5%. 



IV. RESULTS 

In this section we present GFMC spectroscopic overlaps, SFs, and ANCs and compare 
them with those obtained by some other methods and with experimentally deduced values. 
A convenient parametrization of the overlaps is also provided. 

Table [T] shows computed and experimental binding energies for A < 7 nuclear states 
relevant for this work. The errors shown in parentheses are only the Monte Carlo statistical 
errors; in addition there may be systematic errors from the GFMC algorithm of the order of 



1-2% [27|, l38|, |40] . Compared with previous GFMC works, the statistical accuracy on binding 
energies is much better thanks to the high sample counts needed for a statistically reliable 
sampling of overlap tails as mentioned in Sec. IIIIl Most GFMC results were obtained for 
the AV18+IL7 potential. For this potential, GFMC energies in Table [T] differ by less than 
1% from experimental values; in most cases the difference is less than 0.5%. For the ground 
states of ^Li and ^Li, both Type I and II VMC wave functions defined in Sec. Ill Al were 
used. For comparison of GFMC overlaps with those obtained by other methods, we also 
constructed wave functions of s-shell nuclei bound by AV18+UIX whose energies are also 
shown in Table [H 

Besides getting the absolute binding energies right, it is important for consistent overlap 
calculations to work with nuclear states having the correct one-nucleon separation ener- 
gies. In Tables [III and llllt we show GFMC one-nucleon separation energies computed for 
the desired Hamiltonian H from Eq. ([T]) and for the simplified H' from Eq. along with 
experimental values. Given the agreement between GFMC and the experiment, and the 
related discussion in Sec. IIIH we assume the experimental separation energies when deter- 
mining ANCs from GFMC overlaps and account for possible systematic errors inflicted by 
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TABLE I. GFMC binding energies E for A < 7 nuclei for several 2iV+3iV potentials H. Only 
statistical errors are shown on GFMC values. For p-shell nuclei, the VMC wave function type is 
indicated. The energy shown for was obtained in a stochastic sampling of an explicit solution 
of Eq. ([T|). Experimental values are also shown. 



'^Z{J'^,T) 


Type 


H 


E [MeV] 

GFMC 


exp. 


2H(1+,0) 




AV18 


-2.2247(3) 


-2.2246 


v 2 ' 2 / 




AV18+IL7 
AV18+UIX 


-8.47(0) 
-8.47(0) 


-8.48 
-8.48 


3He(r,i) 




AV18+IL7 
AV18+UIX 


-7.72(0) 
-7.72(0) 


-7.72 
-7.72 


4He(0+,0) 




AV18+IL7 


-28.43(0) 


-28.30 




AV18+UIX 


—28.34(1) 


-28.30 


6uc,/'n+ 1 \ 
rle(U , i j 


T 
i 


AT/"1 Q 1 TT 7 

AV io+ilj ( 


— Z\).zU[o ) 


— zy.z 1 


6Li(l+,0) 


I 


AV18+IL7 


-31.93(3) 


-31.99 


II 


AV18+IL7 


-31.88(3) 


-31.99 


6Li(3+,0) 


I 


AV18+IL7 


-29.59(3) 


-29.80 


6Li(0+,l) 


I 


AV18+IL7 


-28.43(3) 


-28.43 


^He(|-,i) 


I 


AV18+IL7 


-28.66(3) 


-28.83 


^Li(i-,1) 


I 

II 


AV18+IL7 
AV18+IL7 


-39.08(3) 
-39.00(3) 


-39.24 
-39.24 


^Li(i-,i) 


I 


AV18+IL7 


-38.88(3) 


-38.76 


^Be(|-,i) 


I 


AV18+IL7 


-37.61(4) 


-37.60 



this choice. 

Tables IIVI and |V] summarize, respectively, spectroscopic factors for s- and p-shell nuclei 
with A < 7. In these tables, the errors on VMC and GFMC SFs are only statistical; possible 
systematic errors on GFMC values were estimated in Sec. Illll to be 2-3% or less. 

In the s-shell, where VMC wave functions are already very good approximations of true 



eigenstates j27|, the Si/2 spectroscopic factors in Table HVl change by less than 2% between 
VMC and GFMC. VMC, however, seems to underestimate the small amount of d-waves in 
A = 3 nuclei by about 15%. Not surprising, s-shell results for AV18+UIX and AV18-I-IL7 
are similar because both 3N forces, IL7 and especially UIX, were tuned to reproduce the 
binding of s-shell nuclei. The experimental information on SFs in the s-shell is scarce; for 
the (^H I ^He) overlap, the GFMC value of 1.61 is in good agreement with an experimental 
value ~1.4-1.6 deduced from electron scattering 

The A = 3 nuclei contain 1.5 S,T = 1,0 nucleon pairs IIqI, Isj. The sums of s- and 
d-wave SFs for overlaps of '^H and ^He with a deuteron are about 1.32 and 1.33, respectively. 



If these values are interpreted as numbers of deuterons [51i, it appears that about 90% of 
T = pairs in A = 3 nuclei are in the deuteron state. We conjecture that the slightly bigger 
(^H I •^He) SF compared with that between and is due to the fact that it is easier to 
compact a spatially extended deuteron into '^He, which is somewhat bigger than the triton. 

Using a similar probabilistic interpretation for ^He, where all one-nucleon (neutron or 
proton) SFs add up to 2 [l| and those for overlaps with A = 3 nuclei in Table [IV] are about 
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TABLE II. GFMC core-valence separation energies in the s-shell for several 2A^+3A^ potentials H 
and for the corresponding H' from Eq. ([5]). The combined statistical errors on the core and the 
parent energies from Table [I] are all < 0.01 MeV. Experimental values are also shown. 



parent 


core 


H 




B [MeV] 










GFMC, H 


GFMC, H' 


exp. 




2H(1+,0) 


AV18+IL7 
AV18+UIX 


6.24 
6.24 


6.24 
6.24 


6.26 
6.26 


3He(r,l) 


2H(1+,0) 


AV18+IL7 
AV18+UIX 


5.49 
5.50 


5.49 
5.50 


5.49 
5.49 






AV18+IL7 


19.96 


19.94 


19.81 


^He(0+,0) 


AV18+UIX 


19.89 


19.88 


19.81 


3He(l^i) 


AV18+IL7 
AV18+UIX 


20.71 
20.63 


20.69 
20.62 


20.58 
20.58 



1.6, we find that about 80% of all nucleoli triples in ^He are in the triton or '^He state. 
Because both A = 3 nuclei are spatially more extended than ^He and the triton is somewhat 
smaller than ^He, we conjecture that it is more likely to find a triton than '^He inside ^He, 
which is reflected by a slight difference in the corresponding SFs in Table IIVI In the s-shell, 
these conjectures are plausible since the s-wave SFs in Table IIVI account for such a big 
fraction of the total spectroscopic strength. 

In Table ilVt the AV18+UIX values labeled HH are those of the hyperspherical harmonics 
method [47]; in particular, the overlaps for A = 3 parent nuclei are from \lm and those for 
He are improved and updated versions of previously published results [15|, H. The HH 



TABLE III. Same as Table jn] for ^ < 7 p-shell nuclei bound by AV18+IL7. (Negative) Positive 
values mean that the A-body nucleus is particle (un)bound relative to the core. The VMC wave 
function type is indicated. 





parent 




core 






B [MeV] 






T) 


Type 




Type 


GFMC, H 


GFMC, H' 


exp. 


^He(i- 




I 


6He(0+,l) 


I 


-0.54(4) 


-0.47(1) 


-0.45 






I 

II 


6He(0+,l) 


I 


9.88(4) 
9.81(4) 


9.82(1) 
9.74(1) 


9.98 
9.98 


^Li(f- 




I 

II 


6Li(l+,0) 


I 

II 


7.15(4) 
7.12(4) 


7.16(1) 
7.06(1) 


7.25 
7.25 




I 

II 


6Li(3+,0) 


I 


9.49(4) 
9.42(4) 


9.45(1) 
9.37(1) 


9.44 
9.44 






I 

II 


6Li(0+,l) 


I 


10.65(4) 
10.57(4) 


10.60(1) 
10.52(1) 


10.81 
10.81 






I 


6Li(l+,0) 


I 


6.95(4) 


7.00(1) 


6.77 








6Li(l+,0) 


I 


5.69(4) 


5.68(2) 


5.61 


'Bed" 




I 


«Li(3+,0) 


I 


8.03(5) 


7.98(2) 


7.79 








«Li(0+,l) 


I 


9.18(4) 


9.13(2) 


9.17 
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TABLE IV. SFs for overlaps between s-shell nuclei for a given 2A^+3A'^ Hamiltonian H and angular 
momentum channel Ij. The statistical errors on VMC and GFMC values are less than 1 in the last 
digit shown; systematic uncertainties in GFMC values are 2-3% or less, as discussed in Sec. Illli 
Also shown are the results of the HH method and experimentally deduced values. 





core 




1 . 

h 


}J 






kJ 












HH 


VMC 


GFMC 


exp. 










AVI 8-l-TT,7 




1 39 


1 30 






'H(l+, 


0) 


Sl/2 


AV18+UIX 


1.30 


1.32 


1.30 




"3/2 


AV18+IL7 




0.0194 


0.0224 










AV18+UIX 


0.0225 


0.0195 


0.0223 










Sl/2 


AV18+IL7 




1.32 


1.31 




3He(l+,l) 


'H(l+, 


0) 


AV18+UIX 


1.31 


1.32 


1.31 




d-i/2 


AV18+IL7 
AV18+UIX 


0.0222 


0.0190 
0.0191 


0.0221 
0.0221 




4He(0+,0) 




h) 


Sl/2 


AV18+IL7 
AV18+UIX 


1.60 


1.64 
1.63 


1.61 
1.61 


~1.4-1.6 [7] 


3He(r 




Sl/2 


AV18+IL7 
AV18+UIX 


1.58 


1.62 
1.62 


1.60 
1.60 





results were converted to our conventions for overlaps. The agreement between GFMC and 
HH is very good, as can be seen from SFs in Table [IV] and from Fig. [6], where the actual 
(^H I ^H) and (^H | ^He) overlaps are shown. Small discrepancies exist between GFMC and 
HH overlaps at short distances, and the (i-wave ANCs in A = 3 nuclei differ; the sources of 
these discrepancies remain unknown. An equally good agreement is obtained for overlaps 
involving ^He. 

As examples of our s-shell overlaps obtained for the AV18+IL7 potential, we show in 
Fig. [7] the (^H | ^H) and (^H | ^He) overlaps. Also shown in these figures are WS fits of 




r [fml r [fm] 



FIG. 6. (Color online) Comparison of the Si/2 and the d^/2 (^H | ^H) (left), and the Si/2 (^H | ^He) 
(right) overlaps obtained for the AV18+UIX potential within VMC, GFMC, and HH. Only statis- 
tical errors are shown on Monte Carlo overlaps. 
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r [fm] r [fm] 



FIG. 7. (Color online) VMC and GFMC Si/2 and 4/2 (^H | ^H) (left) and si/2 (^H | ^He) (right) 
overlaps obtained for the AV18-I-IL7 potential. Only statistical errors are shown. Also shown are 
WS fits of GFMC overlaps. 



GFMC overlaps. Although the absolute overlap signs are rather arbitrary since they depend 
on the detailed construction of our wave functions, the order of vector couplings in Eqs. ([H]) 
and (Q, and the isospin Clebsch-Gordan coefficient in Eq. ([H]), relative signs of different 
overlaps for a given parent-core combination may be important for interference effects. In 
our sign convention, the s- and d-wave overlaps in A = 3 nuclei have opposite signs consistent 
with experimentally deduced negative ratios of their ANCs shown in Table I VII Note the 
slight dips in s-wave overlaps near the origin that are consistent with a depletion of one-body 
densities of s-shell nuclei at short range [27| due to a repulsive potential core. 

VMC and GFMC SFs between A = 6 and 7 p-shell nuclei bound by AV18+IL7 are listed 
in Table |V] along with shell-model predictions and experimentally deduced values. The 
Cohen and Kurath [4^] shell-model values (SM) in Table |V] were corrected for center-of- 
mass motion effects by a factor A/{A — 1) j^, 52] to make them comparable with our results, 
which, as mentioned in Sec. Ill CI are free of such spurious contaminations. In addition, 
a square of the isospin coupling coefficient from Eq. ([8]) was factored into SM values in 
Table |Vl Because it is not clear whether experimental SFs, often compared by their authors 
with original shell-model values of were corrected backward for center-of-mass effects 
or not, we quote them as they appear in the literature. 

In the p-shell, the structure improvement provided by GFMC over VMC is known to be 



more significant compared with the s-shell 27|. This is reflected by SFs in Table IVj that 



change by as much as 15% between VMC and GFMC. The results involving the ground 
states of ^Li and ''Li for which both Type I and Type II VMC wave functions were used 
demonstrate that, by closing the gap between VMC SFs as big as ~20% down to less than 
3%, the GFMC method is rather insensitive to starting trial functions. 



In Table El the experimental SFs from 26| are based on electron scattering data, the 



other experimental values were deduced from hadronic reactions. The experimental SFs 
from jl^] listed in the table are only relative; the authors of that work concluded that their 
data analysis cannot be trusted to provide an absolute determination of cross sections, and 
only relative SFs for several reactions were provided upon a renormalization by a somewhat 
arbitrary factor of 0.32, making them comparable to theoretical VMC values available at the 
time. In ^], SFs deduced from electron- induced proton knockout on ^Li leading to the two 
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TABLE V. SFs for overlaps between A = 6 and 7 p-shell nuclei for the AV18+IL7 potential in 
angular momentum channels Ij . The VMC and GFMC statistical errors are insignificant compared 
with the estimated systematic uncertainties of 2-3% or less for GFMC as discussed in Sec. IIIII SM 
denotes corrected shell-model values (see text). Experimentally deduced values are also shown. 

parent core Ij S 



^ZiJ^'^T) Type ^-^Z(J'^,r) Type SM VMC GFMC exp. 
^He(|~,|) I ^He(0+,1) I P3/2 0.690 0.532 0.565 0.37(7) [21] 

II ^He(0+,1) I P3/2 0.691 If^^ If^^ 0.44(6)1 [23], 0.42(4) [26] 



Pi/2 0.253 0.128 
6He(2+,l) I P3/2 0.212 0.146 



sum 0.466 0.274 0.16(2) [26 



Pi/2 0.338 0.229 0.230 
I I P3/2 0.503 0.480 0.438 
7Li(|",l) 6TY1+ ^-^^^ ^-^^^ 0.74(11)1 [23], 0.73(5) [53] 



Pi/2 0.338 0.211 0.229 
II II P3/2 0.503 0.401 0.428 
sum 0.841 0.612 0.657 0.74(11)1 [23], 0.73(5) [53] 

n ^Li(3+,0) I P3/2 0.646 0.72(14) [5^, 0.58(13) [55J 

II ^^^(""'^) ' ^3/2 0.345 0.19(3)1^ 

Pi/2 0.045 0.069 0.060 

7Li(i",i) I 6Li(l+,0) I P3/2 0.997 0.854 0.759 

sum 1.042 0.923 0.819 1.15 [56], 0.90(9) [53] 

Pi/2 0.338 0.229 0.225 

''Li(l+,0) I P3/2 0.503 0.480 0.438 

7Be(|",i) I sum 0.841 0.709 0.663 

^Li(3+,0) I P3/2 0.646 0.500 0.457 

^Li(0+,1) I P2./2 0.345 0.221 0.210 

^ Values obtained from (ctcxp/cdwea) x 0.32, see text. 



lowest states of ^He were found in perfect agreement with VMC values available at the time; 
in Table W\ the GFMC SFs involving ^He(O^) still agree perfectly with experiment, but the 
agreement for the reaction leading to ^He(2"'") has been spoiled by an error in the VMC 
code discovered after the original work j26| had been published. The (^He(2"'") | ^Li(3/2~)) 
overlap involves an unbound state of ^He for which GFMC does not find a stable energy; 
instead, the method produces the core nucleus with an ever-increasing radius by breaking it 
gradually into "^He and two neutrons. Consequently, the GFMC SFs for this overlap steadily 
decrease and as such are absent in Table |Vl Given the experimental and systematic GFMC 
uncertainties and bearing in mind all the issues related to the meaning of spectroscopic 
overlaps and factors mentioned in Sec. [H we conclude that the GFMC results in Table |V] 
seem to support newer determinations of (relative) SFs. 

We do not find a significant difference between SFs of mirror nuclei ^Li and ^Be; also, SFs 
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FIG. 8. (Color online) VMC and GFMC V\I2 (left) and P3/2 (right) (^Li(l+) | ^Li(3/2-)) overlaps 
obtained for the AV18+IL7 potential. Only statistical errors are shown. Also shown are WS fits 
of GFMC overlaps. Superimposed as dashed lines are the "density" functions (R x r)^. 



between the ground state of ^Li and the J'^, T = 0^, 1 isobaric analogs of ^He and ^Li differ 
just by a factor of 2 due to the difference in isospin coupUng coefficients in Eq. ([H]). Our 
SFs for overlaps between A = 6 and 7 nuclei suggest a substantial quenching of shell-model 
values by as much as 40% (except for the weak pi/2 channel of the (®Li(l"'") | ''Li(l/2~)) 
overlap). 

As examples of p-shell overlaps, Fig. [8] shows the pi/2 and J93/2 (®Li(l+) | ^Li(3/2^)) over- 
laps obtained for the AV18+IL7 interaction with GFMC wave functions originating from 
Type I VMC functions (see Fig. IHfor another example). Also shown in the figure are WS 
fits of GFMC overlaps. The figure illustrates that the overlaps may change from VMC to 
GFMC even when their SFs remain the same (pi/2 channel); on the other hand, sometimes 
the change in SFs is more due to a renormalization than a shape change (P3/2 channel). 
To better appreciate where the changes in SFs come from, we superimpose in Fig. [H] the 
"density" functions [R x r)^ whose integral in Eq. flTT]) gives the SFs. 

The (^He | ''He) overlap in Table IVl is particularly challenging for both theory and experi- 
ment because it involves a parent nucleus that is particle unbound by about 450 keV relative 
to the core. In our calculations, ^He is treated as a bound state and the GFMC propagation 
yields a stable energy and separation energy shown in Table [T] and Table IIIII Despite the 
bound-state approximation to ^He, GFMC substantially improves the overlap tail compared 
with VMC, as can be seen in Fig. [9] where the desired asymptotic form is represented by the 
scattering Coulomb function plotted at 90° phase shift. 

ANCs, extracted from GFMC overlaps by a fitting procedure outlined in Sec. IIIII are listed 
in Tables |Vl] and IVIII for s- and p-shell nuclei along with experimentally derived numbers 
and those of some other realistic methods. Systematic errors on our ANCs were estimated 
in Sec. Illll to be ^5%. As was mentioned earlier in this section, the absolute signs of GFMC 
overlaps and ANCs are not meaningful, but relative signs in different overlap channels for a 
given parent-core combination matter. 

In the s-shell, GFMC and HH overlaps obtained for the AV18-I-UIX force agree well 
(see Fig. [6]), and so do ANCs in Table I VI I except for the weak ci- waves in A = 3 nuclei 
(and consequently for the ratios of d- to s-wave ANCs) for which the two methods differ 
for unknown reasons. The HH ANCs assume the actual separation energies for a given 
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FIG. 9. (Color online) VMC and GFMC P3/2 (''He(0+) | ^He(3/2-)) overlaps obtained for the 
AV18+IL7 potential. Only statistical errors are shown. Also shown is the irregular Coulomb 
function G{r)/r for the experimental separation energy representing the scattering asymptotics at 
90° phase shift. 



Hamiltonian, which are close to the experimental values, and so the ANCs would not change 



dramatically if experimental separation energies were used instead. In [20|, an integral 
method (IM) was used to compute ANCs in A < 9 nuclei from AV18+UIX VMC wave 
functions assuming experimental core- valence separation energies; in Table IVI[ these ANCs 
are labeled VMC-IM. For the dominant s-wave channels, VMC-IM results are in a good 
agreement with ours. On the other hand, the difference between GFMC and VMC-IM 
ANCs for the weaker ci-waves in A = 3 nuclei, larger than 15%, is outside of the systematic 
error bars set on GFMC ANCs; this discrepancy is most likely due to VMC wave functions 



TABLE VI. ANCs for GFMC overlaps between s-shell nuclei for a given 2A^+3iV Hamiltonian H 
and angular momentum channel Ij. Systematic uncertainties on GFMC values are 5% or less as 
discussed in Sec. IIIII Also shown are the results of the HH and VMC-IM methods, and experimen- 
tally deduced values. The errors on VMC-IM are only statistical. For A = 3 nuclei, ratios of d- to 
s-wave ANCs are also shown. 



parent 


core 


h 


H 






C 










HH 


VMC-IM 


GFMC 


exp. 








AV18+IL7 
AV18+UIX 


2.15 


2.13(1) 


2.14 
2.14 


2.11(3) f57],2.07(2) [13|, 1.87(14) [58| 


^H(iM) 


2h(1+,0) 




AV18+IL7 
AV18+UIX 


-0.0925 


-0.0979(9) 


-0.0848 
-0.0842 








Cd/Cs 


AV18+IL7 
AV18+UIX 


-0.0430 


-0.0460(5) 


-0.0396 
-0.0393 


-0.0418(15) [59] 








AV18+IL7 
AV18+UIX 


2.16 


2.14(1) 


2.10 
2.10 


2.10(16) [58], 1.76(11) [60] 


^He(lM) 


2H(1+,0) 


d-i/2 


AV18+IL7 
AV18+UIX 


-0.0865 


-0.0927(10) 


-0.0762 
-0.0794 








Cd/Cs 


AV18+IL7 
AV18+UIX 


-0.0400 


-0.0432(5) 


-0.0363 
-0.0378 


-0.0389(42) [59] 


4He(0+,0) 




31/2 


AV18+IL7 
AV18+UIX 


-6.47 


-6.55(2) 


-6.45 
-6.49 


7.36(19) [58],6.70(50) [58], 5.44(15) [61] 


3He(i + ,i) 


■51/2 


AV18+IL7 
AV18+UIX 


6.36 


6.42(2) 


6.45 
6.49 


6.77(51) [58], 6.52(49) [58] 
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TABLE VII. ANCs for GFMC overlaps between A = 6 and 7 p-shell nuclei for the AV18+IL7 
potential in angular momentum channels Ij. Systematic uncertainties on GFMC values are 5% or 
less as discussed in Sec. IIIII Experimentally deduced values are also shown. 





parent 




core 




h 




C [fm-i/2] 




^ ) 


±ype 




Tvnp 












I 

TT 
11 


6He(0+,l) 


I 


VZ/2 


3.52 

O.DO 


2.48 [62] 






T 
1 


6Li(l+,0) 


T 
1 


Pl/2 
^3/2 


1.73 
9 9Q 

2.87 


1.26-2.82 [63] 






TT 
x± 


TT 

J. J. 


Pl/2 
P3/2 


1.85 
9 90 

2.87 


1.26-2.82 [63] 






I 

II 


6Li(3+,0) 


I 


P3/2 


3.50 
3 fi3 


2.06-3.00 f63] 






I 

II 


6Li(0+,l) 


I 


P-i/2 


—2 3Q 
-2.46 


1.71-2.62 [63] 


^Li(i- 




I 


6Li(l+,0) 


I 


Pl/2 
PZ/2 


-0.57 
-2.85 
2.91 




^Be(r 




I 


6Li(l+,0) 


I 


Pl/2 
P3/2 


-1.70 
-2.20 
2.78 










6Li(3+,0) 


I 


P3/2 


-3.49 










6Li(0+,l) 


I 


P3/2 


-2.58 





employed by VMC-IM that, as shown in Table HVl in terms of SFs, underestimate the strength 
of d-waves in A = 3 nuclei. Our ANCs are in a fair agreement with some experimental values 
but differ from others. Experimentally, particular emphasis was put on the ratio of d- to 
s-wave ANCs in and ^He most precisely inferred from tensor analysing powers; the 
experimental values listed in Table |Vl] agree well with ours, including the relative negative 
phase. 

As shown in Table IVIIl ANCs for overlaps between A = 6 and 7 nuclei are constrained 
rather poorly experimentally. In general, ANCs in the lower p-shell are determined mostly 



from hadronic processes, such as transfer or knockout reactions; see [13|, |20] and references 
therein. In Table IVIIt we show the full range of experimental ANCs from 63| for overlaps 
between the ground state of ^Li and the lowest states in ^Li. Overall, given the experi- 
mental uncertainties, it is hard to compare GFMC to experiment. ANCs calculated within 
VMC-IM for AV18-I-UIX are in a broad agreement with our AV18-I-IL7 results, though with 
some notable differences especially for the P3/2 (^Li(l"'') ] '^Li(3/2~)) overlap. Based on sim- 
ple analytical arguments, the ratio of (^Li(l+) ] ^Be(3/2-)) and (^Li(l+) ] ^Li(3/2-)) ANCs 
involving mirror nuclei ^Be and ''Li was predicted to be 1.02 for both pi/2 and P3/2 chan- 
nels [g^]. We obtain 0.95(7) for the pi/2 channel and 0.98(7) for the P3/2 channel assuming 
5% errors on GFMC ANCs, using averaged Type I and Type II C^Li(l+) ] ^11(3/2")) ANCs 
from Table IVIIj and disregarding the signs of involved ANCs. 
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TABLE VIII. Fitting parameters from Eq. ([THD for GFMC overlaps obtained for the AV18+IL7 
potential in angular momentum channels Ij. All p-shell overlaps originated from Type I VMC 
wave functions. The precise values of Vws are adjusted to reproduce the experimental separation 
energies Bexp- 
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h 


Vws 


Rws 




/3 


P 


Vso 


Rso 




Bexp 


^ZfJ'' T) 






[MeV] 


ffml 


[fml 


0.00 


[fml 


[MeV] 


[fml 


ffml 


[MeV] 




2H(1+, 


0) 


''1/2 


— 172.88 
-2732 90 


0.56 
-1.15 


0.69 
0.91 


1.15 
1.16 


0.64 
0.41 


2.95 


2.15 


0.14 


6.2572 
6.2572 


^He(r,^) 


2h(1+, 


0) 


•5 1/2 
"^3/2 


-179.94 
—8155 10 


0.54 
—2.19 


0.68 
0.91 


1.13 
1.03 


0.64 
0.35 


1.47 


2.07 


0.06 


5.4934 
5.4934 


4He(0+,0) 




\) 


Sl/2 


-202.21 


0.93 


0.66 


0.87 


0.81 








19.814 


3He(i+ 


^\) 


Sl/2 


— zuu.ao 


U.oo 


u.oy 












on c;78 
ZU.D ( o 




6He(0+,l) 


P3/2 


— oo.yo 


z.bo 


U.yo 






i.zb 


U.yi 


n n7 


y.y /DO 




6Li(l+ 


0) 


Pl/2 
P3/2 


-41.80 
-69.55 


3.18 
1.89 


0.85 
1.17 






2.45 
2.13 


1.35 
2.36 


0.17 
0.21 


7.25 
7.25 




6Li(3+ 


0) 


PZ/2 


-62.98 


2.35 


1.18 






1.72 


0.55 


0.09 


9.436 




6Li(0+ 


1) 


P-i/2 


-59.39 


2.64 


0.97 






0.93 


0.91 


0.22 


10.813 




6Li(l+ 


0) 


Vl/2 
P-i/2 


-33.71 
-65.00 


3.39 
2.04 


0.31 
1.15 






1.44 
1.14 


1.09 
0.88 


0.03 
0.03 


6.772 
6.772 


'Be(|-,i) 


6Li(l+ 


0) 


Pl/2 
P-i/2 


-39.45 
-72.22 


3.32 
1.85 


0.76 
1.11 






3.47 
2.62 


1.44 
2.53 


0.39 
0.31 


5.6055 
5.6055 


6Li(3+ 


0) 


Pi/2 


-59.20 


2.52 


1.07 






1.08 


0.79 


0.03 


7.7915 




6Li(0+ 


1) 


Pi/2 


-59.49 


2.64 


0.96 






1.06 


1.05 


0.03 


9.1685 



Finally, we present in Table IVnT] potential parameters appearing in Eq. (fT8|) that provide 
good fits to GFMC overlaps. These fits were used to extract ANCs. They also allow our 



results to be easily used in DWBA or CCBA reaction calculations [45|, |46[. Besides setting 
the core-valence separation energies to their experimental values, no other constraints were 
imposed in the fitting procedure. In the s-shell, we found it necessary to introduce a short- 
range repulsive Gaussian potential term to reproduce dips in GFMC overlaps near the origin 
as seen in Fig. [71 in the p-shell, such a term is not necessary. Good fits to ci-wave overlaps 
in yl = 3 nuclei seem to require very strong repulsive potential cores and negative radii of 
the central WS part. Also, for d- waves the first few radial bins were omitted from the fit 
because they are rather uncertain. In order to reproduce the results presented in this work, 
the WS fits need to be normalized to SFs from Tables IIVI and |Vl Our GFMC overlaps are 
available from LH 



V. SUMMARY AND CONCLUSIONS 

We have reported Green's function Monte Carlo calculations of one-nucleon spectroscopic 
overlaps in nuclei with mass numbers A < 7. The calculations have used wave functions 
derived from a realistic Hamiltonian that reproduces well the low-lying spectra of light 
nuclei. The overlaps are extrapolated from mixed estimates between VMC and GFMC wave 
functions, and they are extended to regions beyond the nuclear surface, where the GFMC 
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method may not be accurate, by means of WS fits to the overlap interior. A goal of this 
work is to provide a consistent set of spectroscopic overlaps, SFs, and ANCs in light nuclei, 
obtained from our currently best GFMC wave functions, that can be used as structure input 
in analyzing existing or future experimental data. 

The comparison of SFs and ANCs with experimentally deduced values is obscured, as 
mentioned in Sec. [H by the model-dependent way these quantities are extracted from exper- 
imental data and by the issues related to the meaning of spectroscopic overlaps. For many 
overlaps, it is hard to judge the agreement or disagreement between theory and experi- 
ment because of the nonexisting or confiicting experimental values. Our (relative) SFs seem 
to broadly support more recent values deduced from hadronic processes, and some agree 
particularly well with values provided by electron-scattering experiments. Our calculations 
reproduce, within error bars, the experimentally well-deduced ratios of d- to s-wave ANCs 
in y4 = 3 nuclei. We observe a rather substantial (up to 40%) quenching of GFMC SFs when 
compared with the traditional shell-model. The GFMC imrpoves the VMC overlaps, but 
the corrections to SFs are sufficiently small indicating that the VMC values used to analyze 



experimental data in |21|-|26[] were reliable. 



The GFMC overlaps presented in this paper are available from [19|]. The overlaps for 



somewhat heavier nuclei up to A < 10 will be the subject of a forthcoming paper. 
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